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ABSTRACT 


Four techniques forthe numerical solution of partial 
differential equations and eigenvalue problems were 
investigated. Typical problems considered were elliptic 


partial differential equations of the form 


De Bes = f(x,y), (1) 


or 


U, + U. + AU = O, (2) 


Mire appropriate boundary conditions are specified so 
that the problem is self-adjoint. 

The four methods are relaxation, Galerkin, Rayleigh- 
Ritz, and dynamic programming combined with Stodola's 
method, for eigenvalue problens. 

mie resules indicated that for eigenvalue problems 
Semaxation or dynamic programming modified 1S TG be 
preferred usually and for partial differential equations 


Galerkin or dynamie programming is preferred. 
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0 ELTON 


Initially for this thesis it was planned to investa 6 
methods for finding eigenfunctions and eigenvalues, with 
@ Particular interest in the oscillation in basins such as 
harbors and bays. The report by Angel [10] introduced 
anew mothod, that of dynamic programming, for the solu- 
ion of some partial differential equations. This method 
meemea promising and it was then extended here to the solu: 
Eon of other partial differential equations. It also 
EE. possible το invert a linear differential operator 
mem nence apply Stodola's method in finding eigenvalues 
en eigenfunctions. This led naturally to a comparison 
With several procedures already known to try to compare 
eemvergence, speed, and accuracy, by applying them to the 
۳۰۰۰۰.٣67 of several rather simple problems. 

It is assumed that the reader has some familiarity with 
meee techniques of replacing differential equations by dif- 
mmence equations and the standard technique of separating 
Variables in linear partial differential equations. The 
regions and their boundaries were assumed to be "nice," 

Nov excessively irregular. Some knowledge of calculus 
of variations also is desirable but not necessary; ref- 
erence [5] provides more than adequate background. 

Four problems were divided into two categories. In the 
first category were eigenvalue problems. Two typical ones 


Were 





Subject 


U... sepa OH 
XX yy 


to the constraint 


0-1-0 


τα Δ. ۲1) 


on OA, (15523) 


where SA is the boundary of the domain A; and second 


In the second category were equations 


equation 


subject 


and the 


subject 


and 


Problems in the first category were numerically solved 





V*U s AU 

to the constraint 
270 

Doy) = - 0 
gin 


NZ + “yy E 


to the constraint 
ο AS I 
biharmonic equation 
V*U = f(x,y) 
bosthereonstrarnes 


U(x,y) = ney) 


=) 
on 


in A, (1:3) 


orm ôA. (1.4) 


Such as Poisson's 


in A, (HS) 
ON OA, πη) 
in A, (1.7) 
on 6A (το) 
on SA. ο 


by a relaxation method, the Rayleigh-Ritz method, and 


dynamic programming combined with Stodola's method. 





Problems in the second category were solved by Galerkin's 
method, the Rayleigh-Ritz method and dynamic programming. 
The domaín used throughout this paper for comparisons 


was the unit square, with the constraint 
UGS) 50 on SA. (1:10) 


Computations were also carried out for L-shaped and trian- 
war regions and for other boundary conditions. 

In the relaxation method and dynamic programming method 
in solving eigenvalue problems the initial estimates to 
the eigenfunctions had a special property. The first ap- 
proximation to Ui was picked so that it had no negative 
value in the domain. The approximation to U2 was picked 
So that it had negative and positive values in the domain 
۲۰٠٠ addition it was orthogonal to Ui. The initial ap- 
proximations to Us was picked in the same way as the one 
Neml except that it had to be orthogonal to both Uj and 
U2. It may be difficult to make a suitable choice for some 


Of the higher modes. 





II. RELAXATION METHOD 


For many years relaxation techniques have been used 
to solve differential equations with and without the 
aid of computers. They are basically iterative proce- 
dures in which a new approximation is obtained from a 


previous approximation and its residuals. 


A.  DERIVATION OF EQUATIONS 

In this section a typical problem is posed and solved 
by a relaxation method. 

Suppose the problem to be solved has the following 


Borm 


Ze. cu. d TN (2.1) 


where the unknown function Z must satisfy the differential 
Son in a simply connected region A in the xy plane, 
ου. ο t>0. In addition the function Z is required to 


vanish at points on the boundary 6A of the region A, 
Z(x,y,t) = 0 on SA. (22) 


A typical problem is that of a vibrating membrane. The 
function Z(x,y,t) denotes the vertical displacement of 
pue membrane. 

Now assume that the displacement has a representation 


er che form 


Wc vov) Wc) Ux). (6.5, 





When Eq. (2.3) is combined with Eq. (2.1) the variables 
may be separated and a new equation is obtained of the 


Form 


" ) 
کے‎ O OEE E Ca) 


Ens is equivalent to two equations 
qur oss و‎ C255) 
and 


ο 
δε t Nr SU, (2.00 


each with a parameter A. Further U must satisfy the 


Boundary condition 
U(x,y) = O on 6A. (2270 


This is a typical eigenvalue and eigenfunction problem: 
“ind the values A, or is and the associated functions 
U ου U, ιτ: pav (ο ο) οπα μις eni E 

ΛΠ... بی ےد ۸22۸ کا کر‎ are Called sine 
eigenvalues. With each eigenvalue is an eigenfunction 
UL: jesin so la 71 
different eigenvalues are orthogonal on the region A. 

Each eigenfunction Ux is also called a mode. 

Consider a thin elastic membrane of a particular form, 
such as a very thin uniform sheet of rubber. Assume that 
pue membrane is made fast at the boundary, while it is 
tightly stretched over the region with uniform tension. 
Also assume that damping is negligible. Then if an interior 


Megion of the membrane is pushed in a direction perpendicular 


10 





Es the plane of equilibrium, 10 becomes αὐθιο νο πο ο 


curved surface. The resulting area can be computed as 
- [[VT* U. 7 * U 7 
S> OS U? dxdy, (2.9) 
Α 


Assume that U, and E. are very small; Eq. (2.8) becomes 


approximately 
S= ff(1* XUL? * XU, ?) dxdy. (2.9) 
A 


The increase in the area of the membrane due to the distor- 


tion is therefore approximated by 


ôS 


{1 + U? + 3U?) dxdy - ff dxdy (2.10) 
A A 


1 ES + Uva? ανα. 
A 


Hence the potential energy of the membrane in the deflec- 


ped position is 


E TS U) dxdy, το ΠΠ; 
Α 


Where VU IS Ορος tensions accumed to be constant over the 
region. 

Now consider any particular eigenfunction or mode. 
It follows from the solution of Eq. (2.4) that the deflec- 
tion is a periodic function of time and may be expressed 


in the form 
Ley toe) A US (2.12) 


except for a phase shift. Thus Eq. (2.11) can be rewritten 


ا 





as 


ΜΠ U?) dxdy sin? At. (215) 
A 


Ehe maximum value of the potential energy is 


PE ax "o /2 [[(U,* * U,*) axáy. (2.14) 
A 


The kinetic energy of an element dm = pdxdy of the membrane 


is 

tpdxdy(U,)? = %pdxdy(U?,? cos?,t), (2.15) 
mirë p denotes the mass per unit area of the membrane. 
Dem ore, the kinetic energy of the vibrating system is 


KE = 425 f [UR dxdy Bas (ARL) 
A 


and the maximum value of the kinetic energy is 


- L2 2 

(ο)‏ 10+ 9م 
Α‏ 

IE". 1S assumed that the energy is constant then, tne 

Maximum values of the potential and kinetic energy are 


Eu for individual modes, and therefore 


5i*p [[U'dxdy - v/2 ff(U, + U,*) dxdy (2.18) 
A A 
or, 
M BS + U) dxdy (2.19) 
2 = 
o ffU?axay 
A 


T2 





The first eigenfunction Ui is the function which minimizes 
tnis quotient and Re e Ve 11% is the corres- 
ponding value for the quotient. The next eigenfunction 
is the function U2 which minimizes the quotient in the 
space of functions orthogonal to U1, and A2* is the corres- 
ponding value of the quotient. The third eigenfunction 
U3 minimizes the quotient in the space of functions ortho- 
memal to U; and U2, etc. The set Uj, Uz2, ... ís unique 
except that if two or more eigenfunctions have the same 
emana lue, they may be replaced by linear combinations 
۶۰۰۰۰۰ ٦90561٤68. It is و"‎ ΠΡ «λος Ax 

If the mode U is known, Eq. (2.19) can be used to obtain 
A^. An interesting fact is that a rather poor approximation 
Besche first mode U1, chosen to satisfy the appropriate 
Pemmaary condition, will yielé@a surprisitingly good estimate 
for Ai^. This result is apparently due to Lord Rayleigh, 


and the quotient is sometimes called Rayleigh's quotient. 


B. COMPUTATIONAL ROUTINE 

MM this section the results of section A will be used 
to develope a method to solve Eq. (2.6). 

πι ο (2.60) 15 expressed as a difference equation 


then it may be written as 


(U ZU) (U 


. 9 ES . e ES 
der) ΤῊΝ πη 
AL 


h? k2 


1,441 * U4,4-1 7 204 50 


IS (2.20) 
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where h and k denote the x and y mesh size respectively. If 


the mesh sizes are equal then Eq. (2.20) can be rewritten as 


+ 
a Den 


233 ' τ 


(U ) DN) 


Me procedure used to solve the 1 EE 
equation was to pick a function U° which satisfied the given 
Poundary conditions as a first approximation to U(x,y). 
ees function need not be a very good approximation to U, 
and in fact step functions were sometimes used. 

From Eq. (2.19) an approximation to A^ was obtained 
by assuming that the approximation U° was the desired func- 
EU. With this first aPproximation to 1* Eq.(2.21) 

EE cused to obtain an improved estimate of U, The form 


of Eq. (2.21) used to obtain the vth approximation was 


v-1 v-1 ν ν 
ee ar. Go 
ËS, EE 
ke (4 = λα.) 


in which the y-lst estimate of λ was used. By alternating 
(5.19) and Eq. (2.22) a close approximation to A^ and 
U were obtained. The solution converged to the smallest 


eigenvalue A, and the associated eigenvector U,. 


C. HIGHER ORDER EIGENVALUES 

in this section the method is extended to find the 
Merger eigenvalues and eigenfunctions. 

To obtain the larger eigenvalues and eigenvectors the 


Same procedure as that in section B was used except that 


14 





additional equations were added to force the elgenfunctions 
to be orthogonal to those already found. Define Ui Don 

i= 1, 2... MLO ط‎ 9٦٢٦6٠٦ ٢٦٦٥٦۰٢٦۶٢٦ 2٥٣٢٦٢۶٢۶۰٢ ٦٤3 ٦٤71 ۰ ne 

with Aj. All higher eigenfunctions must be orthogonal to 
the lower ones obtained.  Orthogonalization was accomplished 


by subtracting out multiples of the lower eigenfunetions 


already obtained. This was effected by expressing U as 


n 


U = U, - MU. (2.23) 
ار ات‎ ger ὃ ὦ 


where 
nn. U, dxdy 
NEG m S oi d cU (2.24) 


J ffu,? axay 
A 


fer cach eigenfunction desaired a different initial 
epproximating funetion was used. The method gave the eigen- 
values and associated eigenfunctions in numerically increas- 
mee order. The method was very simple and effective for 
lower eigenvalues. 

Ene faul of the method is that convergence was poor 
Eecomputation times were large if the number of intervals 


Mme Doth directions was large. 
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III. DYNAMIC PROGRAMMING [8,10] 


The method of dynamic programming has recently been 
lied to the solution of partial differential equations. 
ter). The difference equation may De regarded as leading 
to a system of linear equations with a large number of 
MONO Was. The method effectively reduces the number of 
unknowns involved so that a number of systems are to be 
ved, each one of much lower dimension. In one case, 
li example, using a grid with N+1 intervals in the region 
each direction, N systems each with N unknowns are solved, 
rather than one system with N - squared unknowns. 

seccion A, the method is applied to the solution Of 
Poisson's equation over a rectangle, following the paper of 
Angel [10]. In section B, the method is extended to a re- 
lated eigenvalue problem by combining it with Stodola's 
famed. in section ©, the method is extended to the soiu- 
mom Of the biharmonic equation, by applying dynamic pro- 
παρ twice. Finally, in section D, the method is 
applied to the eigenvalue program involving V'U - A?U, again 


EN cornbrninge the metnod with Stodola'*s method. 
A. DYNAMIC PROGRAMMING FOR PARTIAL DIFFERENTIAL EQUATIONS 
Consider the solution of Poisson's equation 


UE D pU meN (3.1) 


where U = U(x,y) is subject to given boundary conditions 
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Such as 
Ux mc y on ôA. (5 5) 


It may be noted that Eq. (3.1) is the Euler equation 


associated with the variational problem 


I(U) = Min [f(u,? + U? + 2fU) dxdy (3.3) 
U A y 


۰۰۰۰۰۰۲۱۱١ functron U is chosen Irom the ۹ 757۶ء‎ 
EUH first partial derivatives belonging to L^ over A, and 
satisfying the boundary conditions of (3.2) on SA. 

Bev the region A be discretized by choosing ntl and 
mti equally spaced points in the independent variables 
x and y respectively. Then Eq. (3.3) may be rewritten in 


a discrete version with equal intervals as 


n m 
T U = Mi U. 7 a | , 2 
ES U. D AU و‎ RED (3.4) 
1 ل و‎ 
= 2 2 
+ (U, 5 ΣΙ + fjj U, gh |, 


[ . 
where Uo وا‎ Us o)» ID. وا‎ and QU nb are determined 


from the boundary conditions of Eq. (3.2). Now, if all terms 
Ed. (3.4) involving only boundary values were removed, 
while not affecting the solution, a more convenient form 


AE. (3.4) is obtained 


I(U) = Min O S 
U, . ju m ΜΠΕ EN 
2٦ 
m-1 : m-14 
+ A UE 0 77 NE 
A لول1‎ 9 A 2 لو‎ el 





In vector notation Eq. (3.5) may be rewritten in the form 


n 
I(U) = Min 2 (<QU, ,U,> + <rp Us? + Sp (386) 
U 1Ξ] 
R 
t «Up - Up_y> Un - Ux > + < fpo Up?) 
iis, U, = (Us ,)..., U yu This relation now defines 
R Ro R 
m-1 
a symmetric matrix Q, vectors Pps and fa and a scalar 
Sp by 
2 j=j 
Q = ία DÉI where q. = -1 ROSE =] (2 0) 
ول ل و ل‎ 
O otherwise 
` Up y J=1 
rp = در چا‎ where SE Un j = m-1 
O otherwise 
— 2 
fa (2h و و‎ 
2 2 
ern 


The notation as in STR UR. denotes tne scalar product or 


the two corresponding vectors. The matrix Q remains constant 
while rn and S, are functions of the boundary conditions 
enly. 

Now in order to solve Eq. (3.6) a sequence of dynamic 
programming problems was considered. Let 


n 
Min ) (<QU, , U, + ο. U, (3799 


R' e ہ_ لا‎ 7 


ol ıı 





where ے ولا‎ ne is given by πα. 2(0 ٦ 


can be rewritten as 


ER(V) - um (SQUR , R? t η» Up” + Sp eo 
U 
R 
+ SUR - V,Ug - V> + «f " Ug? + <Q Bii? U 41? 
+t... +S, + SS - U): ut NUT ee U 2). 


This may be done since the minimization over Unis, 
E. can be commuted with the minimization over Up. This 
is a common technique of dynamic programming based upon the 
Permciple of optimality, [8]. This allows Eq. (3.9) to be 
rewritten as 


Fa (V) = iin («QU ern US ES Ge 
U R 


F (V) = <Q 0? u uo EE U ۷ (5) 


and U is known from the boundary conditions Eq. (3.2). 
Since FR(V) is quadratic in V Eq. (3.11) may be rewritten 


in the form 


FAC) E «Αη, V> + <b V> + Cp. (5.15; 


IB bywsubstitutjng from Ba. (3.12) into Eq. (3.10) and 


e hot erëntiating viti rëspëct to Un an expression for 


19 





Un is obtained of the form 


St 
EH 


- (I + Q + A R . 00 


El 


R 


Mis is obtained by substituting related Eq. (3.13) into 
meme (3610) and then combining Eq. (3.10) with Eq. (3782), 


the various quantities in Eq. (3.13) are defined by these 


۶9۲60 as 
A, =I- (I+ @+&8,,,)7> (3.14) 
να: (νυ SS (Ga) CHIS) 
Ch = Chyl t SRT E ον; M Geno) 


with initial values determined from Eq. (3.11) as 


m - I, ΠΤ 
By = ΞΕ... (al) 
E MT PQ) Um U t Sr - se ας (3.19) 
The matrix (I + Q + TED, is nonsingpuTar. 1 . Thustit nas 


Hee ES which may be computed beforehand, since it íis de- 
memaent only upon the type of operator. 
Due to the fact that only the values of Un are desired 


the quantities C, need not be calculated. 


R 
The procedure is to calculate the quantities in 


Eq. (3.14) repeatedly until A: and bı are obtained. Then 


20 





f 


U is found trom πα = Uo: which is known 
from the boundary conditions. Next, Eq. (3.13) is repeated- 


ly solved using the stored values of A, and b 


R Rs and the 


last value of Un as V. 

Thus, the problem was solved by n-1 inversions of sym- 
metric matrices of order m-2. While these matrices may be 
Ere there are efficient computer routines available to 
ΤΈΣ ες. Once the inverses are round, they 
NE stored for future use, since they are based only on 
the geometry of the region A. Thus for several problems 
Er tne same region the inverses may be entered into the 
program as data. Also they have the property that if less 


than ntl grid points are required a reduced number of the 


matrices may be used. 


B. DYNAMIC PROGRAMMING FOR EIGENVALUE PROBLEMS 
im the firs, section of this chapter it was seen that 
Ἱ παπάς programming could be used to solve partial differ- 
EMO! equations. In this section, the program is modified 
to solve a related eigenvalue problem by Stodola's method. 
Consider the problem of the vibrating membrane consi- 
HSred in Chapter II. The differential equation for the 


eigenvalues and eigenfunctions is 


= -)2 j 
Um + Se λ20 Inm, (025205 


where U = U(x,y) is subject to the boundary condition 


Ky) Ξε on SA. (3.219 


AJ 





This problem was solved by two related methods. In 
both Eq. (3.20)is regarded as a special case of Eq. (3.1) 


an which 


f(x,y) 9 -A"U(x,y) in A (3.22) 


and 


\ 
O 


E) on 6A (3.23) 


٣۰۰٠-۶ (۰٠0 6 
The difficulty here was that the function f was 
known only on the boundary and A was unknown. The first 
Step, to overcome this obstacle, was to choose a function 
U (x,y)which satisfied the given boundary conditions. Then 
an estimate of A^ was obtained, say eub OT by using Ray- 
leigh's formula. Next a new estimate of U, say Ul(x,y), was 


obtained using dynamie programming to solve the equation 


05 Uyy ΥΕΝ τ 


By repeating this sequence of Rayleigh's formula and dyna- 
mic programming, a good approximation to the minimum eigen- 
value and the associated eigenfunction was obtained. The 
next two eigenvalues and vectors were also obtained DY TEE 
process of forcing the higher eigenfunctions to be orthogo- 
nal to ones already obtained as was done in relaxation. 

Ihe inverse matrices used in the routine were calculated 

in determining the first eigenvalue and eigenfunction; 

they were then stored so that it was not necessary to recal- 


culate them for subsequent eigenvalues. 


Ce 





2w- 5> 090 ئل‎ Merno] 

The second form is more like the usual form of Sto- 
dola's method and hence a brief of Stodola's technique is 
riven first. 

An initial function Vo satisfying the boundary con- 
ditions is selected. It may be considered to be of the 


following form 
Vor ہ٦۴‎ aU: t... (351909 


where a; is not equal to zero. For convenience Kë was 


normalized; the Lo norm of Vo, ||Vo||, is defined as 


MAX |Vo(x,y)1 7 HVolls (3.26) 
A 


and ||Vo|| is set equal to one. Now consider 


Zi -3 a 

L Vo = πο Ui - 2» U2 -... (en 
and 

Me -Ὢ P χι} 21 

L Vo = ES η apr U2 +... κ. ο 


In Eq. (3.28) it can be seen that the relative size of the 

Memponents U,, U,,... are decreasing by a factor (11/12)%f, 

05/11)... respectively. After a few applications of the 
operator, L !, to Vg, the leading term will dominate. 

meus the functions obtained will approximate V,, except 

mer a constant factor, and the ratio of Successive iterates 


fell approach a constant, *“1/1*. 


aS 





3. Second ROUTE 
This was applied as follows. Let Z η 06 0 ٦ 


he equation LZ a = 8 so that 


E js 
V : (3.29) 


This equation was solved by dynamic programming. It was 
Bound convenient to normalize each iteration. Let 


-Z 
m 


V =. , (3.30) 
ρα 


Men the approximation can be made 


ne oe ہت‎ 


Z 
TEAT 
lu -oseguence of functions Vg, Vi,... converges to Ui. 
The process was terminated when successive approx- 


EER OnsS were sufficiently close together, say, 
ES (5.55) 


[meme e is a preassigned small positive number. The result- 
ins function p iS an approximation Ee 


۲۳۰۰16٤ from Eq. (3.31). 


C. DYNAMIC PROGRAMMING FOR A HIGHER-ORDER OPERATOR 

In section A and B of this chapter it was shown how 
παπάς programming could be used to solve Poisson's equa- 
tion and the vibrating membrane program, respectively. 
mie this section, by enotaer modification to une ٣ 


pne Olharmenic equation can be solved. 
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Assume the problem to be solved is of the form 
VD ٣ in A, Case) 


where U = U(x,y) is subject to the constraints 


U = g(x,y) on 6A (3.34) 
πη on 6A. 


Clearly Eq. (3.33) may be rewritten in the form 


ο Gp eae (32057 


where 


Vo USED $(x,y). (3.36) 


Now, since U and 3?U/8n?^ are both known as the boundary, 
(x,y) may be approximated on the boundary. On a rectangle, 


for example, the following relations determine 9 on the 


boundary 
$(0,y) = -P(O,y) + ps (0,y) (5.51, 
d(n,y) = P(A) ur "yy (n,y) 
6(x,0) = -P(x,0) + U (x,0) 
O(x,m) = P(x,m) +U Caine 


XX 


The usual finite difference scheme may be used to approxi- 


mate U and p and thus the relations of Eq. (3.37) may 


De approximated by 


Gan = Oe p  ) 
ο = -P. E SE (3.38) 
Ga QS n? 


eo 





QU t U‏ - کی تا 

E a o Γ 

Nad n?‏ و 
a 1118 ٦‏ 

30 0 n? 

(e MM αι) 

$ = E T 230۳۷ i,m MI 

im i,m ne 

Eet e, =e a) e لاج او .ولا‎ 
and Po 02 τ 0» Po m? Gs m are known fromthe boundary cons 
3 3 3 3 


ditions. Thus (x,y) is now approximated on the boundary. 
First Eq. (3.35) is solved by dynamic programming to obtain 
the function 9 in the region. Then Eq. (3.36) is solved 


by dynamic programming to obtain the desired solution. 


D. DYNAMIC PROGRAMMING FOR HIGHER ORDER EIGENVALUE PROBLEMS 

This method may be extended to the corresponding eigen- 
value and eigenfunction problem much as before. One such 
Meysoical problem is that of a vibrating uniform plate with 
hinged edges. 


Assume that the differential equation has the form 
ου ο M ο... COR 
and the boundary conditions are 


a 
U = —e=0 on SA (3.40) 


Q2 


The technique of sections B and C may be applied in two 


meeps to the solution of the problem. Let © be defined as 


6 = V7U. (3.41) 
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Then the two equations to be solved" are 
vo =e (la) 


and 


y2yp+l on (3.43) 


subject to the boundary conditions, Eq. (3.40). 
In order to start the routine an initial estimate for 


the function U(x,y), say (U, ) is made. The value of 4 


d 
is estimated as was done in section B. Then Eq. (3.42) 

and Eq. (3.43) are solved by dynamie programming to get the 
next approximation (U, j The same criterion for stopping 


3 


is used as that in section B. 
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IV. RAYLEIGH-RITZ METHOD [4,5,6,9] 


The Rayleigh-Ritz method has been used for many years 
obtain approximations to the solution of partial differ- 
ential equations and eigenfunction problems, In this method 
the problem is posed as a minimization problem, say invol- 
meme an integral. Then some linearly independent functions 
feel! Satisfy the boundary conditions are chosen. The solu- 
Dron is m Dane ar ١ب‎ 6110116-6760 AS 
Finally the coefficients in the approximation are chosen 
memes το effect minimization. This leads to an eigenvalue 
pasblem involving symmetric matrices. 

Mae functions chosen may be, for example, polynomials 
DO” degree, or trigonometric functions. Assume that 
homogeneous boundary conditions are given. Let By 7 o, (x,y) 


be n functions which satisfy these and approximate the solu- 


tion U by 
n 
U(x,y) = ea Ch Oy. (4.1) 
A. PARTIAL DIFFERENTIAL EQUATIONS 
in chis seeelon the solution of 
U_+U_ = f(x,y) ΠΗ Ges) 


XX yy 


is again considered. It is the Euler equation associated 


DEED minimizing the integral 
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iz U” * 2fU) dxdy. qe 
A 


ENbstituting from Ea. (4.1) into Eq. (4.27 and intesretilen 


results in a function which may be written 


μι. (4.3) 


for functions of the form (4.1). The minimization of this 
function and an approximate solution to Eq. (1.3) is thus 


obtained by solving 


36. = 0, k= 1, 2,..., m (4.4) 
The effectiveness of the procedure of course depends upon 


the choice of the approximating functions 9 (xy. 


B. EIGENVALUE PROBLEMS 
The method is also applicable to eigenvalue problems. 
Consider again the problem in Chapter II of finding eigen- 


values and eigenfunctions for the equation 


9-7 
E + ο 12U N 2.5) 


EN cot to the boundary condition 
U(x,y) = 0 on 6A. (ντ 


In many problems Er ene os e 1S the Minimum 
of the ratio of two integrals. This fact was shown in 
Chapter II and Eq. (2.19). 

If an approximation to U is chosen as in section A, since 
each function satisfies the linear homogeneous boundary 


Eondition, the sum shown in Eq. (43.1) satisfies it. If 


E) 








this sum is substituted into Eq. (2.19) the resulting values 
define an upper bound for A, for all choices of the C's. 
Further also, the C's may be chosen to give a least upper 
palmada over the subspace spanned Dy 91, Pa, ..., کی‎ 
He that approximation for U is sübstituted into τα ٣ 


the numerator and denominator become quadratic forms in 


(C peta C) - C. The numerator has the form 
) j ons 5) 
a., C.C., = CAC ES 
151 151 -5 20 
where 
ajj = JI C, d >, δ. ανα. (4.6) 
> ΣΟ τα. 


and the denominator has the form 


n n T hë 
b,,0,C, = CU BC (le) 
i=l jal ^" ۵ 
where 
2.07 Y 4.8 
13 7 10 y ( ) 


These define the matrices A and B. The first eigenvalue 
41 Minimizes the quotient in Eq. (2.19) and hence the mini- 
mum value of 


CE NE 


defines the minimum value of the quotient in the subspace 


penned Dy 6,,..., DN 
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To find this is equivalent to minimizing the ον σαν irie 
form 


=L 


۸7۰۶۹ 2 3 (4.10) 
C 


Seet to the constraint that the second quadratic form 
assumes the value one 


Gace or ene (4.11) 


mae problem may be solved in two steps. First find the eigen- 
values and eigenvectors of B. Let the eigenvalues of B 
be m EU 1m aand the associated EE 


vectors Vi. Let T be the matrix 


gc Lg VoM... V) diag(1/mi,..., W/w,). (4.12) 


Then the transformation 
C= TD (4,13) 
meagmces the constraint (4.11) to the form 


ede nde (4.14) 


Bendition (4.10) becomes 


Min DY‏ = ک5 


in DED, (4.15) 
D 
where 


E - TAT, (4.16) 
subject to the constraint (4.14). 
Hence the value of 41° is the smallest eigenvalue of 


the matrix E. 


Sl 








T PE Uma 


Let Di be the associated normalized eigenvector of E. 
Shen the corresponding minimizing function has coeli cierro 


determined from Eq. (4.13) 


C4 s TP D. (4.17) 


luis valuesof i iS an upper bound for the first elgen= 
Esc. the least upper bound in the space of functions 
spanned by %1,..., E In a similar way the second eigen- 
value of E furnishes an upper bound to the second eigenfunc- 
In. etc. 

ars imp le problems, at least. it seems rasy LO CHOOSE 
the $'s so that a good bound for the lowest eigenvalue re- 
sh It is not clear how to make good choices to get 


good estimates for the higher elgenvalues. 
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V. GALERKIN'S METHOD [4,5,6] 


Sometimes the solution to boundary-value problems in- 
ueuving partial differential equations ean be obtained by 
forming an associated functional in the form of a definite 
ο. οσα] which is to be made stationary. For a problem of 
IN type, the Hayleigh-Ritz method is^often an effective 
Procedure for the determination of an approximate solution. 
meme ver, in many instances it is difficult to find this 
51116010081. in such cases, the Galerkin method is often 


effective. 


A, PARTIAL DIFFERENTIAL EQUATIONS 
Consider the linear homogeneous boundary value problem 


of the form 
EU TORY 0531) 


subject to linear homogeneous boundary conditions. The 
ELO! L stands for a linear differential operator, such 
νυ’. 


Suppose that an approximate solution is taken in the 


form 


U, (x,y) - 


x 
Ia 5 
H 


C oy Gy), (Da2) 


where the coefficients Craters H are constants, CODER 
mined, and, aS in the previous chapter, the functions 


on are picked to satisfy the homogeneous boundary‏ وور 


3 





conditions. The coefficients are dependent ὃν επ ο - 
of functions picked and therefore, must be recomputed if 
cala reser mumber of functlensteare later cesen. Mimete 
Bunction Un satisfies the boundary conditions; dto PPP nO 


in general, satisfy Eq. (5.1). Thus there is a residual, 


LU (x,y) - f(x,y) = R (x,y), (5.3) 


which can be viewed as an error or penalty function. Assume 
the solution for U(x,y) ean be expressed by an infinite 
complete series of these linearly independent functions 


in the form 
0) MC QUE (5.1) 
ki KK 


Now U (x,y) in Eq. (5.2) represents a sequence of partial 
sums which approximate U(x,y). Now if the condition is 
imposed that L(U, ) - f be orthogonal to each function 
$. (x,y) on the demain A, the following set of equations 


are obtained (5.5) 


ff (LU) As Doc ۰ ۰ 00۶9 ۹ + 8 7 
Á 


Or by use of Eq. (5.3) (5.6) 


JI, Coy) S(x,y) day “Men foun e کے‎ 
NS Eq. (5.6) is to hold as n + ~ it follows that 


ἘΠῚ Ra - 0 (sae 


ne 
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Let R(x,y) be an arbitrary function satisfying the 
homogeneous boundary conditions. Since ($,) for k =el Soni 
forms a complete set of functions, constants qomo Banne: 
can be found such that 


pec E OX ور‎ (5.8) 
po 
Thus, 
SÍ 1im R.(x,y) n(x,y) dxdy - 0, (5.9) 
À ne 


for any arbitrary n(x,y), so that by the fundamental lemma 
of variational calculus Eq. (5.7) is true "almost every- 
where". Suppose further that L(U, ) + EUR το ٣ 
feds , by the use of Eq. (5.7). 

Galerkin's method requires that the error function 


R,(X,y) be orthogonal to each of the functions ὁ or that 


k? 
meme. 5) and Eq. (5.6) hold. Now by substituting Eq. (5.2) 


into Eq. (5.5) an integral is obtained of the form 


EI 
[T lan - d $, dxdy - 0, (5.10) 


Ihis results in a system of n linear algebraic equations 

in n unknowns Ctt Ce. Furthermore, the syarem Tomin- 
n 

homogeneous unless the function f(x,y) is orthogonal to 


each (x,y). Eq. (5.10) may be rewritten 
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Ι ΓΩ: 


d Urs, 9, dxdy = Pr dxdy, (5.11) 


In matrix notation this becomes 


BCszF, er 
where 

er TS (5.13) 

B = tb, ,} ju رر تک‎ (5.14) 

a A‏ ے2 ا 
T‏ _ 

۷۷ ۱٦ (5-159 
and 

sa a, δ dxdy (Se 

ia = aia dxdy 0867٦) 


The constants C: are obtained by solving the system Eq. 
Ὁ (5.12). 

o l location, a convenient variation. Une way to get 
an approximate solution for the C's is the following. 
Choose n points of the region A. Evaluate the terms in 
the integral of Eq. (5.10) at these points. This yields n 
equations for the unknowns ns This method, called colocas 
tion, is not generally so accurate but it is quicker than 
carrying out tne integrations to define the coefficients 


Ἑ ο ο πο ομα ο. {5-1Τ). 
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Galerkins method is also useful To οὐ: dire t 
solution of variational problems. Assume it is desired to 


minimize a functional of the form 


I(V) » ffF(x,y V,V.,V.) dxdy. (5.18) 
A y 


5 desired to find an extreme value for the functional 
subject to the condition that V(x,y) is prescribed on the 
boundary SA of the domain A. It is known that if the exist- 
ence of an extremizing function U(x,y) is assumed and that 
the function F possesses continuous derivatives of the 
Second order with respect to its arguements, there results 


the condition 


Q2 


A X X σ᾿ 


for an arbitrary function n(x,y) which has piecewise con- 
tinuous derivatives in A and which vanishes on 6A. This 


memaerived by considering a function of the form 
May) = Uy) + enlx,y), (55508 


and differentiating with respect to €. If n(x,y) is other- 
wise arbitrary, then one form of the fundamental lemma of 


variational calculus requires that 


9 9 H 
ےت‎ ν 0ْ سے‎ ۰ (5.21) 


x y y 


Suppose now that n(x,y) is the kth function 0, (x,y) and 


that an orthogonality requirement is imposed as 


si 





ὃ ὃ _ 
Je (Fu x تک تج‎ dxdy = 0, (5.22) 


tork = log. ۴ 


Finally, assume that the function U(x,y) which effects the 
minimization can be represented satisfactorily by a finite 


series 


n 
U = eo) E 
n e k k 


Eq. (5.22) then defines a system of n algebraic equation 
pube solved for the n unknowns Crores Cy Thus Galerkin s 
method is applicable in the solution of variational problems. 
E cr, much of its value lies in the fact that it is not 


necessarily connected with a variational procedure. 


B. EIGENVALUE PROBLEMS 
Suppose that it is desired to solve an eigenvalue prob- 


fm Of the form 

LU = AU(x,y) Inea (5.24) 
and 

U= 0 on SA. (5-25) 


Assume as in the first section that the function U(x,y) 
can be approximated in the form of Eq. (5.23). The problem 


m sOlved by considering equations of the form 


μ(νωρ - AU, | 9, dxdy = 0 (5.26) 
ceu PP VS T 
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Now by substituting for Ua 0676 d a system is ob- 


tained of the form 


n 
- λΥ. C, = 0 2 
(a, ‚Ar, 420; ο ο 
[Don Je- PON us. ٦ 
where 
aj j = Ja) 0, dxdy (5.28) 
Yn > 0 dxdy (5.29) 


wes Nas nontrivial solutions for the C's if, and only if, 
A= ja. ,- Ay, .[= 0 (5:309 
| i,j Yig! , : 


where A is the determinant of the matrix. Thus the values 
of T may be found by solving the characteristic equation 
(5.30). For each eigenvalue à; there corresponds a system 


of equations (5.28) to be solved for the eigenvector 


fei 


Vi marco "V? (5.31) 


BOr] 1 0, 


Now for each eigenvalue A; this system of n homogeneous 
equations may be solved to give the values RS where 
mors coeificient corresponding to the ith eigenvalue. 

Mus the elgenvalues are obtained, vich ΠΕΙΣ 91 70 
tng eigenvectors. The functions >, should be chosen to 
NESE Wmawever Characteristics the solution is felv te 


have. 
39 
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VI. DISCUSSION OE VARIOUS METHODS 


FOR SOLVING EIGENVALUE PROBLEMS 


In Chapters II, III and IV, three methods were developed 
٠۰۰۰۱۲۱110118 ٦6ط‎ first three eigenvalues and eigenfunctions. 
All three gave satisfactory values for the eigenvalues and 
EN functions, and satisfactory times for the rather simple 
problems considered here. These methods were used to obtain 


solutions of the following two problems: 


O e 
Det Ds = ۸> 00007 Com) 


Ex Cot to the constraint 

VI kua on SA: (655 
and 

το κ! DIAM (6.3) 


mr lect to the constraints 


9 21] 


U(x,y) = E 


> 0 O OAs 8.) 
on 


Imus Chapter the computation and numerical results are 
discussed and compared. 

Mia Ot the methods a Set ormruncetions was needed. 
In the Rayleigh-Ritz method these were the basis for the 
9. 1ج0‎ 00100 1168 functions; in the other two methods they 
mere the initial estimates of the functions. The ones usus 


ally chosen were 
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A کش‎ GR COS 
πο E oN Vay cx. (660 
UP - (x - x?*)(y - y?) - y). ae 


Step functions were also used as first estimates in the 
iterative methods; these increased the number of iterations 
required some, but not much, particularly for the higher 
modes. 

Usually the range of the independent variable was di- 
vided into ten equal sub-intervals, in going to a differ- 
ence equation. In the Rayleigh-Ritz method this did not 
yield sufficient accuracy and it was found necessary to go 
to forty intervals. Twenty-five intervals were also used 
in some computations; the intermediate number was chosen 
because of storage requirements in dynamic programming. 

In the iterative procedures some convergence or stop 
criterion was needed, When a relaxation procedure was used 
together with Rayleigh's formula for estimating the eigen- 
value, computation was terminated whenever the eigenvalue 
ama not decrease by at least 0.002. In Stodola's method, 
the routine was terminated whenever the norm of the change 
¡Mt he function U(x,y) was less than 0.002. 

The relaxation method required the least time for this 
simple problem. Most of the time required for dynamic pro- 
gramming was spent in inverting the matrices. Dynamic pro- 
gramming yielded a very accurate approximation to the 


06 1 7۰ 
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For some reason it was necessary to use forty intervals 
to get the desired accuracy with the Rayleigh-Ritz method. 
When ten intervals were used, the eigenvalues of the matrix 
E were significantly too small, apparently due to errors 
in the integration and transformation routines. There were 
at least two other disadvantages of the Rayleigh-Ritz method. 
Perst it was tedious to program. Second, there may be some 
dmi culties in choosing the functions E particularly 
if some of the higher modes are desired. The results of a 
set of computations are given in Computer Output 5 , in- 
i e thë three eigenvalues, the associated coefficients 
nene values of the corresponding functions at various 
Points. The routine is shown in Computer Program 5. The 
necessary matrix transformations and solutions for the 
eigenvalues were carried out using programs TRED2 and TGL2 
respectively, [7]. 

Ihe simple relaxation method of Chapter II had the ad- 
ce of being the simplest to program and to run. It 
had the disadvantages that terminal convergence was slower 
than in the method of dynamic programming and if a large 
number Of points were involved the computing time increased 
ercacly. The results of a set of computations are given 
Mmimcomputer Output 1. The routine is shown in Computer 
٣۲۲۰٠27۶۹۳ 1. 

Dynamic programming had the advantage of yielding very 
accurate values in a small number of iterations. The dis- 


advantages were that it was relatively difficult to program 
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and that it took quite a bit of time to generate the inver- 
ses Of the matrices required. -Computer Output EE c e 
results by this method and the program used is in Computer 
Program 6. 

The values of the eigenvalues obtained by the three 
methods are compared in Table I. In Table I are the eigen- 
values for the first differential equation, the difference 
equation with ten intervals, together with the results of 
the computations, the number in parenthesis by an entry 
indicates the number of intervals used in the computation. 
The eigenvalues for forty intervals are intermediate between 


those for ten and those for the differential equation. 


Differential| Difference Rayleigh Relaxation Dynamic 
Equation Equation(10 


4.44289 4.42463 4.4110 (25)/4.42486(10)| 4.42442(10) 


— | | 4 44751(40)|4.44611(40)| 4.43753(25) 
6.92711 7.23029(40)|6 





Table I. 


Comparison of First Three Eigenvalues for Eq. (1.1) 


The three methods gave close agreement for the first 
eigenvalue but tended to differ on the next two. 

For the differential equation of higher order, only 
relaxation and dynamic programming were compared; these 


comparisons are shown in Table II using ten intervals. 
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Differential Difference Relaxation 
Equation Equation(10) 


۳ ۲92) 
49.348050 


Dynamic 


Programming 


19.599956 19.60767 19.57777 
47.985220 48.00555 47.98473 


49.348040 47.985220 48.02386 47.98474 


able TE; 


Comparison of First Three Eigenvalues for Eq. (1.3) 


Computing times were similar, around twelve seconds 
for each. Computer Output 2 shows the results for the re- 
location method, and the program is Computer Program 2. 
Computer Output 9 shows the results for dynamic programming 


and the program is Computer Program 6. 





VII. SDISCUSSLONOISUEDHE, ΟΙ Τμ ΟΕ ἡ 


PARTIAL DIFFERENTIAL EQUATION BY VARIOUS METHODS 


In Chapters III, IV, and V, three methods were developed 
for solving partial differential equations. These methods 


were used to obtain solutions of the following problems: 


gu + US 2(x* * y X y) d (Tl) 


ο to the eonstraint 

Nat on SA: (7.2) 
and 

vu - 8 in A, SN 
Ect to the constraints 


U(x,y) = 0 on 6A 


2ly -y?) Bor E E (7.4) 





-2(y -y?) fone ο ۷٦ 








2(x -x?) ror y =o 





-2 (x -x°) daje < 


In this chapter the computation and numerical results are 


discussed. 
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In Galerkin's method, an approximation for the function 
O{x,y) satisfying the constraint Eq. (7.2) was made by choos- 
ing suitable functions ®ı (x,y), ®2(x,y) and %3(x,y) to satis- 
the constraint in Eq. (7.2). The approximation for U(x,7) 


was 
UC) OCA I ERA 
ST TOG TCC TLC 


The values of C1, C2 and C3 were obtained by techniques 
described in Chapter V. Computer Output 3 shows the values 
ENEMCOSGS Co, Cas and the function U(x,y) at various points. 

In this method the region A was sub-divided into ten equal 
sub-intervals for each independent variable for the inte- 
gration routine. A second approximation for U(x,y) was made 


for this method as 


Ulx,y) = Cıdı + C202 * C49, το) 


= (x -x?)(y -y?)(xC€i + yl2 + xyCa). 


The purpose of this was as follows. There generally is some 
Skill and art involved in choosing the functions وتوہ"‎ Qa 
DON In fact in the first choice $4 i1S actually the desired 
ndon. The second set of functions was chosen so as to 
get some feel for the consequences of a poor choice of the 

o, 5. The values of C,, C,, C, and the function at various 
points is shown in Computer Output 10. 


The numerical solution of Eq. (7.1) by dynamic program- 


ming with the constraint of (7.2) yielded the values in 
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Computer Output 6. Again the region was divided as was 
ehe in Galerkin's method. 

In the Rayleigh-Ritz method, an approximation for the 
solution, U(x,y), satisfying Eq. (7.1) was made by choosing 
suitable functions %1(x,y), ®2(x,y) and %3(x,y) to satisfy 


the constraint (7.2). The approximation for U(x,y) was 
U(x,y) = Cıdı + C202 + C393 E) 
= (x - x?)(y - y?)(C1 * xyCo * x?y?C4). 


Mie values of C1, C2, and C3 were obtained by techniques 
described in Chapter IV. Computer Output 4 shows the val- 
ues of Cı, C2, C3 and the values obtained for U(x,y) at 
meet Ous points. In this method it was found necessary to 
Enbedivide the region into forty equal sub-intervals for 
mre independent variable in order to get satisfactory 
mecuracy. 

The best approximation to U(x,y) was obtained by Galer- 
kin's method using Eq. (7.5), where % was the desired 
unetion. There was no error and the method was able to 
detect that this was the case. However, when Eq. (7.6) 
Was used the maximum error was eight thousandths. The time 
required to solve the problem by this method was 0.57 
Seconds. 

Ihe problem was Solved by Gymamic pregramming in three 
seconds with accuracy to Six digits. However, over half 
of the computer time was spent obtaining inverses, which 
could have been fed as data from solving Eg. (6.1) in this 


particular case. 
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The Rayleigh-Ritz method required just under ten seconds 
and had a maximum error of two thousandths. 

Because of the time required and the accuracy obtained 
by Rayleigh-Ritz, only Galerkin's method and dynamic pro- 
gramming were used to solve Eq. (7.3). 

Galerkin's method obtained the same values and accuracy 
urhe solution of Eq. (7.3) as it did in the solution to 
Eq. (7.1) and took the same time. 

Dynamic programming obtained the same degree of accuracy 
co t vine Eq. (7.3) as it did in solving Eq. (7.1) and 
took three seconds. The results are shown in Computer 


Output 8 and the program is Computer Program 7. 





VETI ECON EUSTON 


The three methods considered for eigenvalue problems 
yielded satisfactory values for the eigenvalues and the 
eigenfunetions. Generally, the relaxation method seemed 
to be most satisfactory. It was straight-forward to pro- 
gram, and it was faster then Rayleigh-Ritz and dynamic 
۰۰٠٠٠00108. It converged rapidly even if step functions 
Mere msed on several different tests figures, such as the 
L-shaped ana triangular regions. 

The dynamic programming method converged in the same 
Amber of iterations as relaxation, but gave poorer esti- 
mof the second and third eigenvalues. It of course 
mS much more difficult to program and required more com- 
٣۰٠۰٠٠٢٢٠٢ time due to the needed inverses. 

The Rayleigh-Ritz method seemed to have little to re- 
commend it due to the computer time required. It required 
much finer meshing in order to obtain a satisfactory accu- 
racy. The only advantage it had was that no iteration was 
Required. 

Er the three methods considered in the solution of 
Poisson's and the biharmonic equations Galerkin's method 
was the fastest and gave accuracy comparable to the Ray- 
Meaeh—-Ritz method. It was also the simplest of the three 


methods used to program. 
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Dynamic programming gave the best accuracy generally, 
but it required more computer time than Galerkin's method. 
fayleigh-Ritz had the same difficulties asit daid πα 
the eigenvalue problem and was considered of little use. 
While dynamic programming required more time for the 
ution of both eigenvalue problems and elliptic partial 

Gifferential equations, it was very powerful. It only 
required a knowledge of the function U on the boundary. 
It can be extended to irregular regions pals It obtained 
Meme, food accuracy. Much of the time was spent computing 
eeeminverses. If the same points were used in solving 
5٣76۴۹1 different problems, these inverses could be cal- 
culated once and thereafter entered as data, reducing the 


computer time greatly. 
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COMPUTER OUTPUT 1.1 RZELAXATION) 


EIGENFUNCTION-1 


PATH= 8 OMEGA= 4.424858 
X Y U(X, Y) 
0 +. 1 O +. 1 0.09719044 
D. 1 Ο «3 0.2508801 
D o. 1 O +. 2 0. 3094832 
De L O e 7 0.2526 196 
0 +. 1 0 . 9 0.09705645 
0 +. 3 O +. 1 0.2508801 
0 .3 0 +. 3 0.6512997 
E. 3 0 +. 5 Qe 8061690 
O +. 3 0 7ء‎ 0.6586339 
D. 2 0.09 0. 2531579 
ο +. 5 O +. 1 0.3094832 
0 . 5 0 «. 3 0.8061690 
0 +. 5 O +. 2 1.000000 
B. 5 0 e 7 0.8169372 
0 . 5 0.9 Oe 3131154 
0 « 7 0 +. 1 Qe 2526196 
O0 . 7 0 +. 3 0.6586339 
O +. f 0 +. 5 0.8169372 
O . 7 0 +. 7 0.6658412 
Q . 7 0 . 9 025 017 
0 . 5 0 +. 1 0.09705645 
0 +. 9 0 e 3 02531219 
D . 9 0 +. 5 0.3137754 
0 +. 9 0 +. 7 042551157 
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COMPUTER OUTPUT 1.(RELAXATION) 


EI GENFUNCT ۵۳2/2 


PATH=11 


ο ο ο oo 2 2 OO OO OO OO OO OO OO OO OO 


X 


ο - ο. یہ ہبہ ي لد لد‎ VI VI yn MI Mi Y Y UY U ο ra Be YF سم صم‎ 


ΟΜΕΟΑΞ 


ο ο o o ooo o o oo o oO oo یو و‎ OO OO OO O O 


DW or‏ یں یہ ὦ o D‏ ی یہ D‏ - یىی D zl‏ ےم سرں سں - D‏ تسم رں یں ہہ 


6. 921301 


U(X, Y) 
0.1853049 
0. 4886459 
0.686538 
0.4958898 
De 14895596 
0.3032387 
0.8008109 
0. 9981067 
0.8116993 
043099371 

0.001278793 
0.008393560 
0.01233941 
0.008128799 
0.001685064 
-043083511 
-0. 8038315 
=0. 9953710 
-0.8089499 
- 90992 9f 
-0.1921543 
-0.5006261 
056194299 
20450256251 


De 





EIGENFUNCT TON=3 
PATH=11 


CO OO OO OG OO OO OO OO OO OO Oo OO OO ٥ 


A 


COMPUTER OUTPUT 1.(RELAXATION)} 


(Y سم‎ ta Fa ta ps 


ο می‎ n D - - ہا یہ لے لد‎ VD MI nm یا‎ Lu u ثلدا‎ uU 


OMEGA= 


ω OO OO OO OO OO OO OO OO GO OO OO GO بے‎ 


6.927361 


U (X, Y) 
0.1857997 
De 3045443 

0.002898388 
-0. 3070407 
-0.1916554 

0. 4894530 
0.8029368 
0.01102810 
-0. 8017028 
-0.4998162 
0.6086553 
0.9981196 
0.01235584 
04 99533650 
-0.6194320 
0.4950783 
0.8095817 

0.005508117 
-0. 8110784 
-0 5033802 

0.1890559 
0.3086166 
2.925960 8-05 
-0.3116446 


سم MI A‏ رہ MW We D‏ یہ ο‏ ہم لںا یا - WW Ks O yw We ο‏ رہ 


DIS 





COMPUTER OUTPUT 2.(RELAXATION) 


EIGENFUNCT TON=1 


PATH=17 


O O0 OO oo OO OO OO OO oO OO OO OO OO CH 


X 


مم بینم نم VD UM t oo uy UY e j‏ ی ی 


ο 0 0 7 یہ یہ یہ یہ‎ MI 


OMEGA= 


Oo OO OO OO 0Ο نی‎ oO © و ی ی‎ 0 9 00 0 O 


mW -e.e O sl VI DÄ ο ο رہ‎ VD Y ιο رہ‎ Ut t D SH MU t 


19.60767 


UCXs Y) 


0. 1011795 
0. 25714610 
0.3156426 
C. 2600489 
0.1008759 
0.2580318 
0.6597930 
0.8128256 
0.6714125 
0.2608653 
0.3149233 
0.8084638 

1.000000 
0.8271890 
0.3214785 
0. 2583916 
0.6639176 
0.8214287 
0.6778069 
0.2627081 
Ce 1004035 
02579977 
0.3188558 


54 





El GENFUNCT TON=2 
PATH=16 


O O O O O OO OO OO OO OO OO OO OO © © 


A 


COMPUTER OUTPUT 2.(RELAXATION ) 


ہر یم یم ہم مم ὦὢ Wu‏ لن سا یا ین ی یں VI‏ یہ ο 0 0 A y Y Y‏ 


OMEGA= 


O O OO oO ہہ‎ OO OO OO OO © OO OO OO OO GN CO 


y We Oo Y] VI سا‎ re vi یں یہ‎ LA FY D یہ‎ VI D بپم‎ D SH UU t 


48. 00555 


U(X+Y) 


0.1909801 
0. 4904021 
0. 6059132 
0.4979762 
0.1925782 
C. 3037586 
0.7884125 
0.9826701 
0. 8096838 
0.3130895 
0.003429962 
0.01207077 
0.01897281 
0.01404283 
0.004551671 
“Ce 3064085 
τος 1948888 
-0.9861545 
-0.8114642 
-0.3132687 
Qe 4937618 
-0.5022484 
-0.6227177 


NË 





El GENFUNCT ION=3 
PATH=19 


O o o o o o o o o o o o o o oOo o O O O O O O O 


X 


COMPUTER OUTPUT 2.( RELAXATION) 


in Ul wom ο UJ uy +. Fi B oL مم‎ 


Ὁ nD nN نرہ لہ یہ لد‎ MI u 


OME GA= 


O O O O O O OO OO OO OO OO OO OO OO O O 


nm We ON uUi UY ο sl UU Fi 0 - Vi UY - 0 - UU) ri 


48. 02386 


UCX, Y) 


0.1851881 
0. 3020115 


0.009329475 


-0.3049613 
0219235259 
0.4791670 
0. 7890826 
0.02798434 
-07 52122 
-0.5036968 
0.5902148 
0.9785894 
0.03599443 
-0.9824019 
0292271023 
0.4810330 
0.7979080 
0.02484581 
-0.80 38256 
-0.5073704 
0.1851680 
0.3065956 


0.007537059 
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COMPUTER OUTPUT 3.(GALERK IN) 


VALUES OF C(I) ARẸ 


Cl = 1.000000 
02 = 0 
(3 = 0 
X Y UCX, Y) 


0. 0080999922 
0.01889999 
0.02249999 
0.01890001 

0.008100022 
0.01889999 
0.04409999 
0.05249999 
0.04410003 
0.01890005 
0.02249599 
0.05249999 
0.06250000 
0.05250004 
0.02250007 
0.01890001 
0.04410003 
0.05250004 
0.04410006 
0. 01890007 


oO oo oo oo oo oo cuo OG O و‎ O O O 
e e 

DEI We er Fi Fa t‏ بس Cl Un‏ یں 0 لہ ل لہ لہ لہ 

O O ο ο O O O OO OO OO OO OO OO 9 
5 

O Y) VI We ο ہد‎ An UY t ο Nn MI FY VO - MI UD سم‎ 


Sc 





COMPUTER OUTPUT 4.(RAYLEIGH) 


VALUES OF C(I) ARE 


cl = 1.120410 

C2 = -04.5808935 

c3 = 0.3500372 

X Y U(X,Y) 

0 .25 O «22 0.03816108 
O «25 O «20 0.04937190 
O .25 0 5ھ‎ 0.03599290 
O «20 O «25 Oc0 4937190 
O «20 O «20 0. 06231650 
0 .50 0.75 0.04461559 
O .75 0ٗ 5 0.03999 290 
Ο .75 O 50 0.04461559 
0 .75 ο .75 ο ο πο” 


50 





COMPUTER QUTPUT 5.(RAYLEIGH)] 


EIGENFUNCTTON=1 


THE 


σι 
C2 
C3 


O O O O O OOOO 


VALUES OF CCI) 
- 29.89044 
- 0. 003808264 
= 2.580980 
X 
.25 0 
«25 0 
«25 0 
«20 0 
. 50 0 
«20 0 
«15 0 
«(5 0 
5 0 


OM EGA= 


ARE 


.25 
.50 
«5 
025 
«50 
«5 
«25 
.50 
0 15 


E 


4.447512 


UCX,Y) 


1.073552 
1.401157 
1.028184 
1. 431359 
1.868153 
1.370868 
1.073485 
1.401070 
1.028116 





COMPUTER OUTPUT 5.(RAYLEIGHI 


EIGENFUNCTION=2 


TH 


C1 
C2 
C3 


E 


OO OO OO OO OO 


VALUES OF C(I) 


e25 
«25 
«25 
«50 
«50 
«50 
«5 
ο {2 
«0 


9.664003 
112.0784 


O Oo o Oo Oo O oOo O O 


OMEGA= 


ARE 


«25 
«50 
e 15 
25 
e 50 
e (3 
«25 
«20 
«00 


76230292 


U(X,Y) 


0.3401793 
1.766417 
2.309447 

-0.8598440 

0. 6040002 
1.765845 

— 629948 

-0.8604184 

0.3393202 
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COMPUTER OUTPUT 5.(RAYLEIGH) 


EIGENFUNCT ION-3 


THE 


C1 
C2 
C3 


OO تک‎ OO oo ooo O 


VALUES OF C(I) 


tl 


it 


«25 
«25 
«25 
«50 
«50 
«50 
«75 
«5 
e15 


9,656152 
-112.4433 


OO OO OO OO O OH 


OMEGA= 


ARE 


.25 
.50 
«75 
«25 
«50 
«75 
.25 
50 
.75 


7.247070 


Ut Xs Y) 


=1.630179 
-Oe 8650630 
0. 2322842 
-0.8558773 
0.6035087 
1.761141 
0. 2463630 
1.770327 
2.309128 
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oo OO o o o o o o o o o O o o oOo O O O O O O O O O 


COMPUTER OUTPUT 6. ( DYNAMIC) 


Ὁ N OD vi «ὢῷ c y ہبہ ہا یہ یہ یہ‎ Uu uU Y Y uu Y ہم‎ ba pap 


O O OOo oO O0 O00 000 0000005000500 0Ο 


ο SS wm We ο یں ہے‎ UY m 0 SS ی‎ - ON UW زہ 0 عم نبا‎ MI ο دم‎ 


U(X,Y)] 


0. 008099854 


0. 01889967 
0.02249958 
0. 01889964 


0.008099858 


0. 018895656 
O. 04409914 
0.05249893 
0. 04409913 
0.01889965 
0.02249958 
0.05249896 
0.06249879 
0. C5249900 
0.02249960 
0012089965 
0.04409916 
0.05249896 
0.04409920 
0.01889972 


0. (08099854 


0.01889966 
0.02249962 
02078892072 


0.008099906 


02 





COMPUTER OUTPUT T.ODYNAMICI 


FIGENFUACTION=1 


PATH= 4 


ο ο ο ο OO OO OO OO O OO OO OO ο ο OO یه‎ ο ο ο OO ooo 


X 


ta ga ta 


OO nn ANA تا ی ل لہ ہبہ یہ‎ MI UN uy U uy ο DD UU میم‎ HM 


OMEGA 


O © و‎ ©0 OO OO oO OO OO © oo OO OO OO OO 


6 
wm UJ F2 O بی لد‎ UU rn aO مہ ہم ںں یں لہ‎ wa MI Y tt D SH ی‎ LA Fi 


4.423473 


UCX, Y) 


0, 095542338 


Ce 2500933 
Ce 3091003 
0. 2500929 


0.09554338 


0.2500930 
0.6546414 
0.8090985 
0.6546412 
04. 2500929 
Ce 3091002 
C. 8690987 

1.000000 
0.8090994 
Ce 3691004 
0.2500930 
0.6546419 
0. 8090992 
0.6546427 
0.2500940 


0.09554315 


0. 2500930 
0.3091C10 
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COMPUTER OUTPUT 7.(DYNAMICI 


EI GENFUNCT ION-2 


PATH= 6 


oO ex "Er ο er e e ον e oocco"o'"o"uoococoococ' oooZso نث‎ 


X 


(Ὁ Ώω 4 - رہ‎ nnn Ui in UN U NM Y DD QQ YY UY یم فصل منم‎ o ps 


OMEGA= 


ο ο ο OO ο ο ο ο ο ο ο ο ο ο نی و و و یی مو تی‎ 0 


wm WwW نم‎ OO - An UU η” O AY ی‎ UU 0 یا پیہ‎ A Fé ο Y ہیں‎ ο ta 


6.927172 
U(X,Y) 


0. 1913275 
0.5004815 
0.6183150 
0.5004799 
0.1913269 
0.3092604 
0.8089828 
0.9994567 
0.8089805 
0.3092589 

1 . 95 901 2 * -06 
3.750642'-06 
2.407420'-06 
-4,.548319'-07 
-7.790408'-07 
- (, 3092556 
-0 709089775 
-00.9994535 
-0. 8089827 
-0 e 3092611 
ος ۶۶ +۶ 
-0.5004784 
-N. 6183146 
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EIGENFUNCT 1011-3 
PATH=11 


o o o o Oo o o o o OO OO OO OO OO OO O O 


X 


COMPUTER OUTPUT 7.(DYNAMIC) 


ο ο D ل‎ zl لد ل لد‎ VI uU Ul Ul Ui UJ U Ww U لرا‎ rn یم عم‎ psp 


OMEGA- 


ο ο. ο ο ο ο ο 2 ο ο. ہی‎ ο OO ο oo oo O O 


Υ 


OO SH Ul Uu سم‎ D sel ری‎ LAN kä 0 SS ہیں‎ UU μυ O SY MI LA rm 


dc 


WwW w 


6.927173 
U (X,Y) 


041913257 

C. 3092576 
1۰ے‎ -6 
-0 3092595 
-0. 1913273 
Ce 5604185 
0808981 

ω 06 
270280509915 
-0. 5004808 
0.6183135 
0.9994553 
1:4 8089-06 
= 0. 9994555 
0۔‎ +7 
0. 5004803 

C. 8089837 

4. 195991906 
-0.8089792 
-0. 5004812 
019192659 

0. 3092602 
2622319 °'706 
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ο ο OO OOo ο می‎ ο Oo OO OO OO oO OO OO OO OO 


CƏMPUTER αυτο ο του νο [([ 


kod ہم‎ = ps po ps 


D 0 40 vi 4 - c - بہ لہ‎ UM nn uy U A y uu U 


ae 798955 01 6.2.0 جح یچ‎ O O 6G ο 0 0 ο تن تج نی نف‎ 


w =- ο SS یىی‎ Were ہہ ف‎ MN t ON WW -” ON WwW 5 


M‏ یہ کی 


UCX, Y) 


0.008099731 


0.01889934 
0.02249917 
0.01889930 


0.008099727 


0.01889932 
0.04409828 
0.05249788 
0.04409826 
0.01889930 
0.02249917 
0. 05249792 
0.06249751 
0.05249795 
0.02249918 
0. 01889931 
0.04409831 
0.05249793 
0.04409834 
0. 01889938 


0.008099716 


001829932 
0.02249920 
0. 01889938 


0.008099772 


06 





COMPUTER OUTPUT 9.( DYNAMIC) 


EIGENFUNCTION=1 


PATH= 


O 0 O OO OO OO OO OO OO OO OO OO OO OO 


3 


Ὁ Nn یہ یہ یہ یہ‎ Y Y تا‎ un VI A yn للا‎ Y uy UD UU t یم‎ Fa Pp 


OMEGA- 


O O00 0 OO OO ο ο 0000050000000 


ο UY mnm D y Ul GUY Pr iO Qu Wer SH MI We ο un Y سم‎ 


19537639 


U CX, Y) 


0. 09549332 


0.2500034 
C. 3090197 
0.2500030 


0.09549332 


0. 25060032 
0.6545131 
0.8090189 
0.6545127 
0. 2500032 
Qe 3090197 
0.8090194 

1.000000 
Ce 8090203 
0. 3090200 
0. 2500033 
0.6545139 
0.8090203 
0. 6545147 
0. 2500044 


0.09549326 


0. 2500034 
0 «3090206 
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COMPUTER OUTPUT 9. (DYNAMIC) 


EIGENFUNCTION=2 


PATH= 4 OMEGA= 47 «98485 
X Y U(X,Y) 
O +. 1 0 +. 1 C.1911072 
O +. 1 ο .—3 0.5001712 
O +. 1 0 +. 5 0.6181207 
O +. 1 0 e 7 C. 5001690 
O +. 1 0 .9 0.1911063 
0.3 O cil Qe 3091234 
0 +. 3 0 « 3 0.8090394 
B.23 0 +. 5 0999282773 
0 . 3 GQ aft 0. 8090362 
0 e 3 O . 9 0.3091207 
0 « 2 O «4 1 -9. 384 042*-07 
Ce 5 0 . 3 8.968276" =06 
D . 5 @ . 5 1.008 165'*-05 
U . 5 O . f 3. 054505*-06 
Dn. 5 De H -4.620023'-06 
Wo. 7 O 1ء‎ =C. 3091 261 
ο +. f O e 3 -0. 8690 2187 
077.7 0 55 2049995166 
O e 7 0 . 7 -C. 8690367 
0.7 De H 304690913253 
0 +. 9 0 +. 1 "ος ο LIS 
O +. 5 0 . 3 -0.5001757 
0 . 9 0 . 5 -0.6181255 
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COMPUTER 


GIGENFUNCT ION=3 


PATH= 7 


O O0 00000000 000000000000 


X 


W F^ e e ta HM 


0 OO sell ` لد لد‎ el جا لہ‎ VI Uu An UY wo οὉ y 


OMEGA= 


ο ο SS ο uv oO ο ο ο OO OO ο ο OO نین ین‎ OO OO OO o 


Y 


wn We ο ο سی ی‎ ka iQ gU UD - ο sl ی‎ We 0 -ϱ WwW نم‎ 


OUTPUT 9. (DYNAMICI 


41.9886 


-9, 


-2 o 


De 


3. 


U(X,Y) 


0.1911126 
0. 3091245 
O5BTTIO'-OT 
-0.3091258 
-0.1911120 
Ce 5001729 
0. 8090208 
200 - 76 
-0. 8090366 
-0.5001741 
0.6181221 
0.9998204 
343255'-06 
2009998225 
-0 6181232 
Ce 5001751 
Qe 58090382 
5 [۷ 2:88۲ -6 
048090329 
20 7509 TT 37 
041911140 
03091279 
614680'-06 
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COMPUTER OUTPUT 10. (GALERKIN) 


VALUES OF C(I) ARE 


Cl 
C2 
C3 


145396 322 
1596392 
πο Y 


X Y U( X,Y) 


0.0023756 32 
0.01059511 
0.01862749 
0.02069907 
0.01103618 
0.01059507 
0.03192532 
0.04658196 
0. 04633235 
0. 02294400 
0.01862741 
0.04658184 
0.05916639 
0.05281764 
0.02397245 
0.02069896 
0. 04633217 
0.05281758 
0.042401 38 
0.01732974 


oo o ooo oo oo oo o.o'"oo"o © O O O 
e e 

Hi € ua‏ نم SH zl uU Ul uU oU MI A Y UJ uy WwW t‏ إل لد ل 

CÓ ο ου o. 0O'"oo oo o oo ooo oo o 
e 

O Wwe on Ww tr ON ہیا‎ We ON WW mn 


TU 





COMPUTER ۹۹۶ LS 


RELAXATION METHOD APPLIED TO 


VS) Ui 


mi 

> 

| 

ھا 

9 

27 

2 = 
2 4C 
CÓ et — 

N 
eI $e 
GC: ee e. 

OQ نت‎ — — 
az E D ? 
Gil e E — 
9 C^ Uu uU. ra ο =- 
ا سے‎ i. KM) Loss eec سے‎ | 
WS E ہے‎ = de ett, ~ E 
LACH x DD MY er t —— o = 
COLE T ke AE d ہے ہے‎ O e ——— eum SE ου 
کت‎ OO Pr Cn zu. وسر ےرہ‎ | > 
0 N میم‎ bei ba e ` NH oe NY بے‎ — = ο 
= MY) )سام‎ ei — LINN => 7 C 
OT E Om pri “7 سے سے بف ۶ اسم )20 و‎ umm — O 
ZO oc nam KU De a ےم‎ ٣4 ب صا‎ > y oO 
«r4 "Y y MS ος” ene سے ے |1 ]ا‎ ~ a «(SN 
A = 2 ۰ی کی‎ O "P MC: — 99 9». eo Cit: a N LI 
ο C a ZU eebe? KEE کے سے‎ — — ^e 
12 LU UNO UI ee ae mimm € (N — Al 
Dal Ou οσα Ce gh ou e E rik ee vl — ا‎ om CC? — n 
—HY. ز× مد‎ E Dn, Cl pa Ê ا‎ e <> (Y MN” Cies E Mm 
t e [x OB جچ‎ 200 CARE aw e t) 17,00 ο. — — | on 
2a LO tini pu AD Ze nn ا‎ COLO زم اسم ہے‎ a4 SN) ہم‎ 
ει GC ποιον, AN Ori ج سے‎ 0١ر‎ GORDO ΖΩ -- 0 vis ۰ 
E CG AUS. p Lit) “A σα. on mor H oodd t ee 
OO = ους οι μασ σου E ee m ZN Priti ہے زے ہے ہے‎ x cue ~x کے‎ 
ett ena >zz 2700 E p (MÁS (OD (GN n0 NI Ke NDO uw D c کس ہےر‎ Eu جه صرح‎ 
ly VER. «pet > >OZF ος) Sg eo oem ο >. E- b 0 luc NH tn 
«fex C 05b- © e OTIC LIN « e ee Dre zu X ee ES ! 
ری ص‎ cy p 1:00)» 0X XU Uv tY YF oS mi jt etri لااتات‎ pb As A reap bt CIDE DC UO m~m oe 
O ZS IS EA LAS A OO DHO wrt ee C E ا ہے‎ kerer 11 ہے‎ Zen) E] 
ح2 ہے‎ E (a EE ja S orri Zu DOK e ez TC y e > <7. واس‎ N c. zx wë Z eet ai UU ege p^ CIA ^ nt 
Oo j “ee لہ‎ Ce ου SS Sao ee) E ک‎ te MM ole ا‎ HW HH (De tea"? 
ue لہ ل‎ GEN DO CO NN ey aia = IL ven o ett nia 
I- LL «C > e es TO Ὁ Qeni L0 Ca ll A A IE I E AD DININ 
rale Na e Weil eD Me Tg C CO CC cc Mom wom ye TH ki لم‎ ee y. m 
Cx οἱ NOOO uy Sa I "m Mu zt. )ام‎ ¡Ef eo ro (hill)... .. .. .. io ri ONOK سم‎ 
zb ux zzzeuvee em evui tt ica ο πας ی‎ MTD eT ter ee dE CC) Hr? <î 
ee ر6٥ بہت انت 3ا .تا‎ O6 TEZ ZT DOS I „termin! CT Ni wel nm OI LIA JO LL. C) 
wein usa ا‎ aen || a || oo oot aiai Fig ana 
i CO ee A E »»(χ come hhh e SË Qi” ri. ez ti CC OC OGO ee iL u 
MO c cux Cz IPA T TELE GO Kap H. UL LLiLit U gies e 


1 


un; 
TEE > Tome 
ي06‎ 
20۲۰۰۱٢: "٣٢ 
ze J:z^Ó UNTIL PI 09 
B GIN 
3 ص04‎ 72 
E ADA IRA AU A IN Kee 
ee a (Oe ieee t= 1007, =50609,.03 
SNC IG O; 
ΤΟ» 
Gestalt GT 2° AND Ke: 
IDE sek e $ ei, ei E 
SNE 7:=) UNTIL Nl1-1 DR 
mes =N2(T) UNTIL N3(1)-1 DI 
IN 
ο τι + (C1 211th J 20155) (7٢+1 :)22ي(‎ H 
mae (2 (io Stl) —Z( 1453) )/42)2 -2)7 HS 
EB on 2) Hi?H2; 
OME 22=0M2 3 
ο ae / K-35 M2 =SORT( IM?) 
fee CM A ?$ (N»47) πη ^| 
EIT 
MPTE (INAT CONV Ρατ): 
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COMPUTER PROGRAM 6. 
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